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Abstract 

A family of Zeta functions built as Dirichlet series over the Riemann zeros are shown to have 
meromorphic extensions in the whole complex plane, for which numerous analytical features 
(the polar structure, plus countably many special values) are explicitly displayed. 

Fonctions Zeta pour les zeros de Riemann 

Abstract 

Certaines fonctions Zeta definies sur les zeros de Riemann, par une famine de series de Dirich- 
let, se prolongent a tout le plan complexe en des fonctions meromorphes, dont de nombreuses 
caracteristiques peuvent etre explicitees (la structure polaire, mais aussi une infinite de valeurs 
speciales) . 
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1 Introduction 



This work proposes to investigate certain meromorphic functions defined by Dirichlet 
series over the nontrivial zeros {p} of the Riemann zeta function ((s), and to thoroughly 
compile their explicit analytical features. If the Riemann zeros are listed in pairs as usual, 

{p — | ± ir fc } fc=1 2 v .., {Re T fc } positive and non- decreasing, (1) 

then the Dirichlet series to be mainly studied read as 

oo 

Z(<r,v) = J2(n 2 + v)-°, Rea>±, v > ~n 2 , (2) 

k=l 

extended to meromorphic functions of a G C, and parametrized by v — with emphasis 
on two cases, v = \ and especially v = 0. Their analysis will simultaneously yield some 
results for the variant series 

oo 

3(cr, a) = ^(T fe + a)~ 2a , Re a > \ (and, e.g., | arga| < tt/2). (3) 
fc=i 

Those Zeta functions are "secondary": arising from the nontrivial zeros of a classic 
zeta function (here, C( s ))j an d "generalized": they admit an auxiliary shift parameter just 

like the Hurwitz zeta function (((s,a) = f X]o°( n "'" a )~ s )- Such "£-Zeta" functions have 
occasionally appeared in the literature, but mostly through particular cases or under very 
specific aspects. On the other hand, their abundance of general explicit properties seems 
to have been largely ignored, although it can be revealed by quite elementary means 
(as we will do). And with regard to the Riemann zeros, reputed to be highly elusive 
quantities, those properties constitute additional explicit information: this is enough to 
motivate a more comprehensive treatment (and bibliography) of the subject. The present 
work just aims to do that, in a self-contained and very concrete way, as a kind of "All 
you ever wanted to know about (-Zeta. functions. . . " handbook (without prejudice to the 
usefulness of any single result by itself). 

We begin by developing the background and our motivations in greater detail. 

First, if a Selberg zeta function is used in place of ((s) from the start (assuming the 
simplest setting of a compact hyperbolic surface X here), then the {r^ 2 + j} become the 
eigenvalues of the (positive) Laplacian on X, and the Zeta functions (0) turn into bona fide 
spectral (Minakshisundaram-Pleijel) zeta functions, for which numerous explicit results 
have indeed been displayed (with the help of Selberg trace formulae: cf. (23 for v > \, 
[Ml for v = \, EE! for v = 0)- 

Some transposition of those results to the Riemann case can then be expected, in view 
of the formal analogies between the trace formulae for Selberg zeta functions on the one 
hand, and the Weil "explicit formula" for £(s) on the other hand |rjj|. Indeed, a few 
symmetric functions over the Riemann zeros that resemble spectral functions have been 
well described, mainly V{t) = J2 P ept (OB HE], E3J chap.II). Zeta functions like © have 
also been considered, but almost solely to establish their meromorphic continuation to the 
whole a-plane — apart from the earliest occurrence we found: a mention by Guinand 
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(see also [8J of the series t^ s (= Z(s/2,v = 0)) on one side of a functional relation 
(eq. (|79)l below) arising as an instance of a generalized Poisson summation formula. Later, 
Delsarte introduced that function again (as </>(s) in ^2]) to describe its poles qualitatively, 
displaying (only) its principal polar part at s — 1, as (2n)~ l /(s — l) 2 ; Kurokawa [24 made 
the same study at v — |, not only for £(s) but also for Dedekind zeta functions and Selberg 
zeta functions for PSL 2 (Z) [or congruence subgroups] (then, Zeta functions like (J2J) occur 
within the parabolic components); and Matiyasevich [2Z| discussed the special values 
9 n = 2 2(n, |) (n G N*). Extensions in the style of the Lerch zeta function have also been 
studied (0S|, ES| chap. VI). 

Independently, Deninger [T3] and Schroter-Soule [33] considered a different Hurwitz- 
type family (we keep their factor (2tt) s just to avoid multiple notations), 

e(,,x) d = f (27rr^(x-p)- s (Hea>l), (4) 
p 

mainly to evaluate d s ^(s,x) s=0 (as eq. (jl(Jl|) below); earlier, Matsuoka J2H], then Lehmer 
[26] had focused upon the sums 



2?n = J> _? \ n e N* [= (27r)- n £(n, 1), for n ^ 1]. (5) 



del' 

P 

Here we easily recover £(s, a;) from the other Zeta function (J3J) (but not the reverse), as 



| + y) = (27T)* 3( 1 S) iy) + e -iW2 3( i a> . (6) 



Tie present work proposes a broader, and unified, description for all those (-Zeta 
functions, with a wealth of explicit results comparable to usual spectral zeta functions 

m- 

Tools for the purpose could also be borrowed from spectral theory (trace formulae, 
etc.), but the objects under scrutiny are more singular here (the Zeta functions for the 
Riemann zeros manifest double poles, vs simple poles in the Selberg case); this then 
imposes so many adaptations upon the classic procedures that a self-contained treatment 
of the Riemann case alone is actually simpler. Even then, we cannot get maximally 
explicit outputs for all cases at once (e.g., Weil's "explicit formula" diverges for v < j), 
and our analysis has to develop gradually. 

So, we begin (Sec. 2) by setting up a minimal abstract framework, sufficient to handle 
Z((T, v) (with a permanent distinction between properties in the half-planes {Re a < 1} 
and {Re o > |} respectively). We next obtain a first batch of explicit results for the case 
v = \ (Sec. 3), then for general values of v (Sec. 4). Specializing to the case v = in Sec. 5, 
we reach an almost explicit meromorphic continuation formula for Z(<J, 0) into the half- 
plane {Re a < |}, which immediately implies many more properties of that function, 
and is generalizable to L-series and other number-theoretical zeta functions. In Sec. 6 
we exploit the latter results to sharpen the descriptions of both Hurwitz-type functions 
Z(cr,v) and 3(a,a). Sec. 7 provides a summary of the results; essentially, a Table collects 
the main formulae for v = and |, also referring to the main text like an index. (Which 
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text can in some sense be viewed, and simply used, as a set of "notes" for this Table !) 
Finally, Appendices A and B treat some subsidiary issues: a meromorphic continuation 
method for the Mellin transforms of Sec. 2, and certain numerical aspects. 

For convenience, we recall the needed elementary results and notations [Tj 19*]: 

B n : Bernoulli numbers (Bq = 1, B\ = — |, B 2 = |, . . . ; B 2m +i = for m = 1, 2, . . .) 
(B n (x) : Bernoulli polynomials) 
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Euler's constant: Stieltjes constants : 7 n == hm < 

M-s-oo I ^— ' 



!f (logm) n (logM) n+1 
m n + 1 

m=l 



(7 = 7 « 0.5772156649, 7i « -0.0728158455, 7 2 « -0.0096903632, . . .) (7) 
-En : Euler numbers (E = 1, E 2 = —1, E4 — 5, . . . ; -E2m+i = for m = 0, 1, . . .). 

[i? n is also a standard symbol for quantized energy levels, a concept often invoked purely 
rhetorically about the Riemann zeros; let us then insist that our (present) work is totally 
decoupled from such interpretations, whereas it sees the Euler numbers as truly essential !] 

dcf °° 

Concerning the Riemann zeta function £(s) = Yl k~ s EHl HH Ej we will need its 

k=l 

special values, 

C(-n) = (-1)"-^ (n = 0,l,...); C(2m) = g^yyi^l (m = l,2,...); (8) 

C(0) = -| and C , (0) = -|log27T, (9) 
and its functional equation in the form 

E(s) d ^ f M = C(1 ~ g) F( S ) d ^ f — (10) 

[) F(s)-F(l-sY [S) s(s-l)T(s/2Y [W) 

where H(s) is an entire function, which is even under the symmetry s <-> (1 — s) (this 
expresses the functional equation), and only keeps the nontrivial zeros of C( s )- 

In order to analyze C( s ); we wm be forced to invoke a particular Dirichlet L-series as 
well (associated with the Dirichlet character \a [U EES ELI] ) 

OO 

= f ^(-l) fe+1 (2A;-l)- s [=L( X4 , S )] ; (11) 
fe=i 

/?(s) extends to an entire function, having the special values 

rn / /o\2m+l 

0(-™) = ^p /?(2m+l)= |£ 2 J (m = 0,l,...); (12) 

13(0) = I and /? / (0) = -|log2-logvr + 21ogr(i), (13) 
and a functional equation expressible as 



dcf ^4y«+i)/2 r ^a + l 
where the function S X4 is entire and only keeps the nontrivial zeros of /3(s 



= 5^(1 -a) (14) 



2 General Delta and Zeta functions of order < 1 



2.1 Admissible sequences and Delta functions 

Throughout this work, a numerical sequence {x k } (systematically labeled by positive 
integers k — 1, 2, . . .) will be called admissible of order fi < 1 if: 

(i) < x\ < x 2 < • • ■ , Xk T +°° i 

(or: complex x fc — > oo with | argx fc | sufficiently bounded [SUI2D1, to provide an uncondi- 
tionally valid framework "in case of need"); 

(ii) ^2 \xk\~ 1 < oo, making the following Weierstrass product converge, 

k 

A(z\{x k }) = l[(l + z/x k ) (VzGC); (15) 
k 

it then defines an entire "Delta" function A{z) (we omit the argument {x k } except when 
an ambiguity may result); 

(iii) for z — *> oo, logA(-2) = o(|,s:| W)+5 ) V5 > 0, and it admits a complete uniform 
asymptotic expansion in some sector {| argz| < 9}, of a form governed by some strictly 
decreasing sequence of real exponents {// n }, as 

oo 

logA(z) ~ ^(a^logz + ajj**" (^^oo) (16) 

n=0 

with fi > fii> • • • , fi n I — oo, and < /i < 1 

( "generalized Stirling expansion" , by extension from the = k [22 J; such a uniform 

expansion is repeatedly different iable in z. 
Then, the Dirichlet series 

oo 

Z(o | {xk}) = S ^ x k' 7 (convergent for Re a > fio) (17) 

k=l 

defines the Zeta function of {x^}, holomorphic in that half-plane. The point er = 1 has 
to lie in the latter by assumption (ii), which imposes fj,Q < 1; then (log A) can moreover 
be expanded term by term in eq.()15j) to yield the Taylor series 

OO y / \ 

log AU) = ^(-l)" 1 - 1 ^^^ (convergent for \z\ < Xl ). (18) 
/ — ' m 

m=l 

Motivations: the idea here is to assume certain properties for an entire function A{z) 
of order /i < 1, so as to generate a function Z(a) meromorphic in all of C with poles 
at the /i n , their maximum order r being 2 here (as dictated by the specific form of 
eq.(Jinj)). Such A(z) are very special instances of zeta-regularized infinite products, for 
which much more general frameworks exist (e.g., [22 1221 I2H] )■ However, singularities 
and essential complications definitely increase each time r or the integer part can 
become larger (chiefly, the formalism leaves ([fio] + 1) meaningful integration constants 



5 



undetermined). Efficiency then commands to minimize both latter parameters (subject 
to r > 1 and /io > 0); specially, //q < 1 is simpler to handle than any /i > 1. In this 
respect, spectral zeta functions frequently have /i > 1 (e.g., for Laplacians on compact 
Riemannian manifolds, /io = |x [dimension]) but only simple poles, hence the simplifying 
assumption r = 1 is suitable for them jH7[|M]; by contrast, the present functions Z(a,v) 
will accept the lower value /io = § but will need r = 2. Another difference is that 
for eigenvalue spectra, a "partition function" Ylk e ~ tXk * s a natural starting point; in the 
Riemann case, that type of function (V(t) = Y2 P ept ) exhibits a more remote and contrived 
structure [TU], while Delta- type functions are easily accessible (by simple alterations of 
eq. ffTUjO . Thus, Riemann zeros and standard eigenvalue spectra have several mutually 
singular features, making their unified handling rather cumbersome. 



2.2 Meromorphic continuation of Zeta functions 

The bounds logA(^) = 0(z x ) for z — > 0, and 0(z fx °) for z — > +oo, with /i < 1, imply 
that the following Mellin transform of log A(z), 

r+oo 

1(a) = \ogA{z)z- a ~ 1 dz (/i <Rea<l) (19) 



converges to a holomorphic function of a in the stated vertical strip. 

Then, by standard arguments (see App.A, and EUEZUHl), I( a ) extends to a mero- 
morphic function on either side of that strip: 

a) by virtue of the expansion (jlfi)) for z — > oo, 1(a) extends to all Re a < 1, with 

(at most) double poles at a = fi n , polar parts = M " + - — — — -; (20) 

(a - fi n ) 2 (a - /J, n ) 

b) just because log A(z) has a Taylor series at z — 0, 1(a) extends to all Re a > /x , with 

(at most) simple poles at a = m 6 N*, of residues — (log A)^ m \0)/m\ . (21) 

Now, for /j, < Re a < 1, the integral (J19)) can also be done term by term after inserting 
the expansion ([To]) , giving 



a sin na 



n^)=J2 x k a = z ^)- ( 22 ) 



7T 

fc=l 



The meromorphic extension of 1(a) then entails that of Z(a) to the whole cr-plane as well, 
i.e., to the half-planes {Re a < 1} by a), resp. {Re a > /io} by 6) independently, so that: 
a') any non-integer pole (/i„ ^ Z) of 1(a) generates at most a double pole for Z(a), with 



[i n Sill 7T/i n _ 



/i n sm7r/i n fsmirfi 



Z(/i n +eJ = a^ n — H a M n+ h/i n cos 7r/i n a^ n - +regular part ; 

(23) 



1 
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a") any negative integer pole (//„ = —m) of /(cr) generates at most a simple pole for Z(a) 
through partial cancellation with the zeros of (sin7rcr), with 



(m e N*) : Z(-m + e) = (-1 



i 1 11 



TYX CL ffYi . ^ 

+ (a_ m — mo. 



\m+l 



0(e) 



i.e., residue: Res a= _ m Z (a) = (— l) m+ ma_ m ; (24) 
finite part: FP CT= _ m Z(cr) = (— l) m (a_ m — mo_ m ) (25) 

(we call eqs.([25|) "trace identities" by extension from spectral theory, specially when a_ m = 
0, in which case explicit finite values for the Z(—m) result); 

a'") any pole at a = of 1(a) gets fully cancelled by the double zero of (asm Tier), giving 

Z(e) =a + a e + 0(e 2 ) Z(0) = a , Z'(0) = a . (26) 

V) each positive integer pole of 1(a) (a = m E W, from eq.(j2U) gets cancelled by a zero 
of (sin7TO"), giving 

Z(m) = (-l) m - 1 (logA) (m) (°)/(" ,l - 1 ) ! ( 27 ) 
but this last output just duplicates with the previous Taylor formula (JTHJ). 

Ultimately, all the poles of Z(a) lie in a single decreasing sequence {fi n } n< z^, and have 
maximum order r = 2 under the specific assumption (jlfij) . 



3 Delta function based at s = 0, and the Zeta function 



We now apply the previous framework to the Riemann zeros upon just slight changes 
with regard to the usual factorization of ((s) around s = 0. 

3.1 Basic facts and notations 

The most convenient starting point is the entire function (jlU|) . which keeps precisely the 
non-trivial zeros {p} of ((s), and has the familiar Hadamard product representation ^3] 

~(s) = e Bs J](l - s/p) e*/', B d ^ [log ~]'(0) = log 2^ - 1 - \ 1 . (28) 
p 

Zeros are henceforth grouped in pairs as {p = ^±irk}k=i,2,... according to eq.((H). Their 



corresponding counting function, J\f(T) = f cardjrfc | Re < T}, follows a well-known 
estimate 77 (T), as [HIIH 



jV(T) ~ 77 (T) = ^- log f - - 1 (T — > +oo). (29) 



2tt 



2tt 



Thus, the sequence of zeros itself could be admissible of order 1 at best; fortunately, 
a transformed sequence {xk} and its Zeta function will immediately arise: 



- fl +r fc 2 , and Z(a) d ^ Z(a \ {x k }) = (30) 



l 



A' 



(the latter series will converge for Re a > ~, again by the estimate 

Indeed, once the zeros have been reordered in pairs, it first follows that !%.\ = f ^2 p _1 = 

p 

J^Xfc -1 = Z(l) (convergent sums), and then, H(s) = e [ jB+ ^^ 1 ^ s [1 + s(s — l)/xf.]; the 
k k 
parity property H(s) = 2(1 — s) thereupon imposes Z{1) = —B, hence 

Z(l) = & = -[log S]'(0) = - log 2v^F + 1 + \ 7, (31) 

a classic result ([H], ch.12; [T4"| . Sec. 3. 8). All in all, the product formulae and functional 
equation boil down to 

((S).a|r n _ C(l-g) fcY \ def 718/2 1 ^ 

= 1 = FO^) i (S) = a (a-l)r(V2) J ' (32) 

where A(A | {x k }) = E(s) with A d = s(s - 1) : (33) 

i.e., S(s) has been rewritten as an infinite product A(A | {xk}) which is manifestly even 
(under s < — > (1 — s)), and will qualify as a Delta function of order fiQ = \ in the variable 
z = A, a much simpler situation than /i = 1 (naively suggested by eq.(|2BJ))- We now 
derive the ensuing properties of Z(u) (as compiled in Sec. 7, Table 1). 



3.2 Properties of Z(a) for Re o < 1 



As basic initial result, the sequence {x^ = f | + r^ 2 } is admissible of order ~: it obviously 
fulfills assumptions (i)-(ii); the function S(s) is entire of order 1 in s, hence | in A; finally, 
a large-A expansion of the form (J16)) for logA(A | {xk}) is easily obtained as follows. In 
eq. ((32) ? for s — > oo (|args| < vr/2), — logF(s) can be replaced by its complete Stirling 
expansion and logC(s) = 0(s~°°) can be deleted, giving 



lo 7T S 1 S 1 

logA(A) ~logA-^-s+^(logs-log2-l)--log| + -log27r+J] 



r> o2m- 1 

JD 2m^ 2m+l. 



2 2 & ^2m(2m-l) 

(34) 

whereupon s is to be substituted by the relevant solution branch of s(s — 1) = A, namely 



s - 2 



\ + \A+| ~ A 1/2 + § + V -f^l£M\W-». (35) 

n=l v ' ' 



The resulting A-expansion then has all the required properties, with the exponents {fi n = 
^(1 — n)} — giving the order /x = \ — and leading coefficients 

5i/2 = I oi/a = -|(! + log27r) (36) 

a = | a = | log87r (37) 

a-i/2 = ^ a-i/2 = -^log2vr - etc. (38) 

This expansion can be computed to any power /i n in principle; still, a reduced general 
formula for a Mn looks inaccessible this way; by contrast, all the a Mn with fi n ^ arise from 
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the single substitution of eq. ((3*3)1 into the prefactor | of log s (~ | log A) in eq.("2U), giving 

ai/2- m = —r— ■ 5-i-m = 0, m = 0,1,2,.... (39) 

m! T(— m + 3/2) 

Now, by Sec.2.2a'-a'"), a large-A expansion for log A(A) translates into explicit prop- 
erties of Z(a) for Re a < 1. Here, Z(ct) gets a double pole at each .half-integer \ —m (m e 
N), with principal polar term i— ^— (| — m)ai/ 2 - m / (c + m — |) 2 , and is regular elsewhere; 
the leading pole <J = \ has 

r „ i 11 l0g27T 1 

full polar part = — t- r . (40) 

P P 87r(a-|) 2 4tt a-i V ; 

At a = 0, eqs.(jHZJ) and ("26|) deliver two explicit values, 

Z(0) = 5o = g, Z'(0) = ao = | log87r (» 0.806042857); (41) 

the latter makes the Stirling constant (the value exp — [Z'(0)], or regularized determinant) 
also explicitly known for this sequence {xk}. 

Further quantities, tied to the yet unspecified general coefficients a^ n , will also acquire 
explicit closed forms: the polar terms of order l/(cr + m — |) and the finite values Z{—m) 
for all m € N, see Table 1; but we will need a more indirect approach (Sec. 6.1). 

3.3 Properties of Z(a) for Re a > \ 

Now by Sec. 2. 26'), Z(a) is holomorphic in the half-plane Re a > |, and the values Z(n) 
for n = 1, 2, ... lie in the Taylor series of log A(A) at A = 0, which can be specified here 
through eq.(|3*2*jl. 

We first expand log A(A) in powers of s, where A = s(s — 1). Eq.("2l"]) directly implies 

oo w 

logA(A) = logH(s) = -E (J?n = ^ P "" ); (42) 

n=l p 

then, from eq. ("3~2~|) . \og[( (s) / F(s)] and log[C(l — s)/F(l — s)] (expressing the functional 
equation) respectively yield the two Taylor series 

- logv^F s + logr(l + s/2) + log(l - s) + log[-2C(s)] 

= - io g ^ + [-1. + 1 tW s »] -±ls + ± 02£« s ., (43) 

. 2 ^-^ 2 n n J z — ' n ^— ' n! 

n=2 n=l n=l 

log(l - s) + \^(s - 1) + logr(^)] + log[-<(l - s)) 

n=2 n=l K ' 
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The 7^ in the last line are cumulants for the Stieltjes constants 7 n of eq.(J7J) [T] (cf. also 
the rj n in ref. jSHl, Sec. 4), i.e., 

n=l V '' n=l ^ 

(7g = 7 = 7, 7? = 7i +| 7 2 « 0.260362078, 7^ = 7 2 +2 7 7i +| 7 3 « 0.034459088, . . .). 

The identification of the three series (pH^ - lpFfj) at each order s n now yields a countable 
sequence of 3-term identities: the first one just restores the result 3F\ = —B as in eq.(JHH); 
then, the subsequent ones likewise express the higher 3f n in two ways, 

2^ = 1- (-1)"2— C(n) - (l °f KI)( 1 t (Q) = 1 " (1 " 2-")CH + T-^rry 7^, n = 2, 3, . . . 

(n — iy (n — iy 

(46) 

That short argument subsumes several earlier results. The rightmost and center expres- 
sions in eq.(|46|) amount to formulae for 3f n by Matsuoka .29, and Lehmer (|26j. eq.(12)) 
respectively; the implied relations between the derivatives C (0) an d the Stieltjes con- 
stants 7 n were also discussed in [2J E] , together with Euler-Maclaurin formulae for the 
£(™)(0) which parallel the specification of the 7 n in eq.(J7j). 

As for the values Z{m) themselves, they are given by eq. (j2?|) now using A = s(s — 1) 
as expansion variable, i.e., 



\m—l 



Z(m) 



(m - 1)! 



1 



2s - 1 ds 



logS(s) (m = l,2,...). (47) 

- s=0 or 1 



But alternatively, Z{m) = ^ m + [a finite linear combination of the {<^} n =i,...,m-i]j an d 
vice-versa: as shown below, 

0<e<n/2 V 7 1=0 V 7 

(Z(l) = 3Ti, Z(2) = ^2 + 2^1, . . .). The jT n being already known from eqs.jSD and 
Z(m) then reduces to an explicit affine combination (over the rationals) of B = [logS]'(0), 
C(n), and either (log |C|) (0) or ^n-i for 1 < n < m (see Table 1), as stated earlier 
by Matiyasevich [27|. (More recently, Z(l) and Z{2) also got revived in studies of the 
distribution of primes |3lH l2~5]). 

Since ref. j^Z] uses eq. ()48|) for Z(m) without mentioning any proof, we sketch one. First, 

if x = f p(l—p), the expansion of the identity log[(l— s/p)(l — s/ (1— p))] = log[l— s(l— s)/x] 
in powers of s yields [p~ n + (1 — p)~ n ]/n = ^2 0<i<n / 2 (— l) e ( n ~^ e ) x~ n+e / (n — £) for n = 

m 

1,2,.. .. By recursion, this must invert in the form x~ m = Yl V m ^ n [p~ n + (1 — p)~ n ]; 

n=l 

m 

then, ^2 V m ^ n p~ n has to be the singular part in the Laurent series of x~ m around p = 0, 

71=1 

resulting in V m;n = fi m ~™~ l \ Now, summing both sets of identities over the Riemann 
zeros {p} yields the stated decompositions (|48j). (We stress that their finite character is 
specific to the Z{m) as opposed to all other values Z(m, v), v ^ j.) 

10 



Note added in proof. - For completeness, we quote two other sets of identities for 
the sums Xi @U|: 




for each odd k > 1 



(a countable sequence of "sum rules" , easy but unreported, which allow to eliminate any 
finite subset of odd values); and the connection to Li's coefficients (cf. [39], thm 2), 

An= f Y P [1 — (1 — Vp) n ] [ w hich allow to recast the Riemann Hypothesis as A n > (Vn)]: 

An = X>W^ ^ = E(-l) i+1 ("V, («=1.2,...). 

3=1 U/ j=l U/ 



4 Generalized Zeta functions i£(cr, i>) = Z(a \ {t^ + v}) 

We begin to discuss the Hurwitz-like generalizations of the preceding case obtained by 
shifting the squared parameters r^ 2 . The obviously allowed translations (x k t— > x k + 
v", v" > —Xi in the real case) preserve the notion of admissible sequences (with their 
values of /io, r); that validates the earlier definition (J2J), as Z(u, v) = Z(o \ {r^ 2 + v}). 

The corresponding transformation of Delta functions (as Hadamard products of order 
< 1) only involves an explicit constant denominator to preserve their specific normaliza- 
tion A(0) = 1, as 

A(A | {r fc 2 + v}) = A(A + v - v' | {r fc 2 + v'})/ A(v - v' \ {r k 2 + v'}). (49) 

Alternatively, we might have opted to normalize Delta functions as zeta-regularized infi- 
nite products, i.e., 

A Z r(A) = exp[-d a Z(a | {r k 2 + A})] CT=0 , (50) 

which are fully translation-covariant, but at the same time less explicit. The overall 
benefit of this normalization is then dubious within the restricted scope of this work; 
but here, it explains a dichotomy between algebraic and transcendental properties of Zeta 
functions, which roughly follows our overall division between {Re a < 1} and {Re a > /io} 
properties, but not quite. 

Covariance implies that if we translate A i— > (X+v ), the expansion of log A zr (A) around 
the invariant point A = +oo can be recomputed to any order by straight substitution, 
yielding explicit polynomials in v as coefficients. When /i < 1 as here, then A zr (A) = 
e -ao ^(^) j[|7| , hence the previous statement holds for the expansion minus its term 
of order A ; i.e., all the shifted coefficients a^iv), a Mn (f) will be polynomial excepting 
ao(v). For logA(A | {r^ 2 + v }) specifically, eqs.tjnSHSHI) v = \) imply that 

«i/2, 0,1/2, oo stay constant (as well as a_i = a_2 = • • • = 0); (51) 

Ol/2-m(v) = ~ ( r(3 (^ /9 , ^7- (52) 

r(-m + o/2)4m! 
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For functions like Z(a,v), the consequences are that their polar parts and "trace iden- 
tities" will depend polynomially on v; furthermore, a single (fixed-v) large-A expansion, 
such as eq.flHU) for logA(A | {r k 2 + i}), suffices to express those v- dependences in full. 
Precisely here, by eq. fl5Tj) : Z(cr,v) keeps its rightmost (a = |) full polar part constant 
(and given by eq. (j4TJ]0 . as well as its value at (Z(0,v) = |); all its other poles (of 
order 2, except at v — 0) keep fixed locations. (As a by-product, any difference func- 
tion [Z(u, v) — Z(cr, vo)} is holomorphic for Re o > — \). We can specify such polynomial 
formulae in the half-plane Re a < | still further, but only later by a different path (Sec. 6). 

By contrast, all formulae for Z(cr, v) in the half-plane {Re a > |} refer to Taylor 
coefficients of log A Z r(f + A) around v finite, which evolve tmnscendentally with v. here 
they will only express in terms of log |C(s)| (or logS(s), from eq. ljlOjO and its derivatives 
at s = -i?; 1 / 2 . The first of those coefficients, log A ZI (v) itself (= — a (t>)), actually yields 
a special value lying at a = 0, by eq.(J5UJ): 

d a Z(<r, v) a=0 = a (v) = a (|) - log A(v - J | {r fe 2 + f }) 

= ilogSvr-logHdit; 1 / 2 ) (53) 

(a result which fully matches eq. (jl(Jl|) below for the Hurwitz-type function £ jTSJ El]). 
Then, the Taylor coefficients of order n > 1 (identical for log A(A) and logA zr (A)) yield 

Z(n, v) = L-l^ — logH(I ± v^) (n = 1, 2, . . .), (54) 
e.g., Z{l,v) = ±\v-V\\og~y{s = \±v^) (v ^ 0), Z (l, 0) = i(logH)"(|) 

(but we are in lack of more reduced closed forms for general (n,v)). 

In summary, the polar parts of Z((J, v) and the special values {Z(— n, v )} n eN ("trace 
identities") have polynomial expressions in v ; whereas {Z(n, v) {n ^ 0) plus 9 CT Z(0, v) at 
n = 0} are also special values, but only computable transcendentally. (This conclusion is 
moreover fully obeyed for typical spectral zeta functions.) 

5 Delta function based at s = |, and the Zeta function 

Z{ct) = Z(a, v = 0) 

An interesting option is now to shift the parameter v from its initial value t in Z(cr), 
to the most symmetrical value f = 0. By eq. (}4T?j) . the Hadamard product (|3^)) becomes 
based at s = | as 

S(i + t) = S(|)A(t 2 |{r fe 2 }) with t = S -|, t 2 = A + i, (55) 

and S(i) = -vr" 1/4 r(|)C(|) (^ 0.994241556) (56) 

(but very little is known about [2H] and we cannot make this constant factor any 
more explicit, contrary to the special case v — j where that factor was S(0) = 1). 
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The factorized representation (}3*2"j) then transforms to 

A(t 2 1 {r, 2 })^| = C(| + 1), D(t) = f C(|) r(|)^V r(f + |). (57) 

(This Delta function is closest to the determinant of Riemann zeros used by Berry-Keating 
for other purposes jl].) 

We accordingly switch to the Zeta function of the sequence {t^ 2 } |Hl[T2], or in short, 

oo 

Zip) = Z(a, 0) = J2 r k - 2a , Re a > §. (58) 
k=i 

Numerically, this new function looks almost indistinguishable from Z(cr) (see App.B). 
(Also, by Sec. 4, [Z — Z){o) extends holomorphically to Re a > — |.) By contrast, the 
meromorphic continuation of Z{<j) will prove to be distinctly simpler than that of Z(a), 
thanks to specially explicit representation formulae in the half-plane {Re a < §}. To 
obtain these, we now switch to a more powerful, special to v = 0, approach (whereas the 
earlier considerations would still describe Z(u), but just to the same extent as Z(a)). 



5.1 The shifted spectrum of trivial zeros 

The factor D(t) in eq. ()57|) has the structure of a spectral determinant built over the 
"spectrum" of trivial zeros in the variable (—t), namely {- + 2k} (D(£) is not exactly 
the zeta-regularized determinant, but again this will not matter here). That spectral 
interpretation can be extended to the factor (1 — 2t) _1 , by treating the pole t = | (of 
C(| + 1)) as a "ghost eigenvalue" of multiplicity (—1). A major role of the spectrum of 
trivial zeros is to make log[D(t)/(l — 2t)] asymptotically cancel logA(£ 2 | {^fc 2 }) to all 
orders when t — » oo in | argt| < tt/2, given that \og((^ + t) decreases exponentially there. 

We therefore expect the spectral zeta function of the trivial zeros (of £(| — t j) to play 
an important role; this "shadow zeta function of C( s )" (f° r short) involves both ((s) itself 
and the partner function /3(s), in the combination 



oo 

z(«) = Ed + 2k ^ s = r [\ (C 1 " 2 "K( S ) + ^) ~ 1 



(59) 



k=i 



(= 2 s ((s, |) in terms of the Hurwitz zeta function). Z(s) has a single simple pole at 
s — 1, of residue », and admits the special values 



Z(-n) 
Z(n) 



n + 1 
(-l) n 2 



B n +i(j) ~ 2 



-n _ 1 
2 



'l_2-«)^+i + 2-"- 1 J E„ 
n + 1 



(n)/5> 



Z(0) = - 




jyp [log r] 



(2 2 ™-l)|5|B 2m |+2^( 2m ) 



- 2 



2m 



2m+l 



l)C(2m+l) + fgJ|E 2m 
and Z'(0) = -|log2- ilog7r + logr(i 



2~ n , n = 0,l,..(60) 



n = 2m 



n = 2, 3,... (61) 



m = 1,2, 



22m+l ^ = 2^71 + 1 



(62) 
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Remark: literally, the framework of Sec. 2 excludes the sequence of trivial zeros (of 
linear growth, and order /i = 1), but the truly relevant function here will be Z(2er), as 
Zeta function of the modified sequence {(| + 2A;) 2 }, which is admissible of order | again. 



5.2 Meromorphic continuation formulae for Z(p) 

We start from a slight variant of the representation ()22j) for Z(cr), obtained through an 
integration by parts upon the Mellin formula (fTTH) (where z = t 2 , by eq. ()55)) ): 



smna def 
Z(a) = J(a), J(a) = 



7T 



(t 2 )- CT dlogA(t 2 |{r fc 2 }) a<Rea<l). (63) 



We next introduce a (regularized) resolvent trace for the spectrum of trivial zeros, 



R(t) 



def 



dt 



logD(t) = \ logvr 



r' 



+ 



(64) 



which has a simple pole of residue +1 at each trivial zero of £(f + 1); a corresponding 
function for the pole ( "ghost" ) at t — | is 



R g (t) = -l/(t - \) (with residue (-1)). 
Then, upon insertion of the factorization formula (fBTjl . eq.(JBHJ) yields 



J(a) 



,-2(7 



-R(t) - R g (t) + ^ + t) 



dt. 



(65) 



(66) 



Now a crucial feature of the case v = is that this integral is also a Mellin transform with 
respect to the argument appearing in the factorized form of ((s) (namely the variable t, 
in eq. (p77]) ). As a consequence, the contribution to J from R(t) (and also R g (t)) can be 
neatly extracted and evaluated, in closed and interpretable form. Because the factor in 
brackets in eq. lJBTIj) is 0{t) at t — (due to the functional equation), J (a) can be split as 



r+oo 

J {a) = J (a) + J T (a), 3(a) = / r 2<T [R(0) - R(t)]dt 

Jo 



(67) 



def 



t 



-2a 



R g (0)-R g (t) + ^ + t)-^) 



dt, 



(this splitting preserves the convergence strip < Re a < 1} for both resulting integrals). 

J r (cr) can be split still further, once its integration path has been rotated by a small 
angle: either +e or — e, in order to bypass the poles of R g (t) and of ^-(| + t) at t — |, 



J r (cr) = Jf(<7) + Jf(<7), 



+ e ±1£ oo 



r^[it; g (o)-it; 6 (0]dt 



(6J 



def 



+ e ±ie oo 



-2cr 



C',i 



^(5+*)-^(i) 



dt. 
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Now, 3(a), Jg(a) can be straightforwardly transformed into Hankel contour integrals 
and then computed in closed form (by the residue calculus), giving 



Z(2a) , , 7r2 2<T e^ ic 



W = -^T-> -fW = — y > (69) 

sin 2na 8 sm 27r<7 

both of which are explicit functions, meromorphic in the whole plane; chiefly, J brings in 
the shadow zeta function (|59*j). 

Then, again upon back-and-forth integrations by parts, J^(a) continue to 

/■+e ±le oo f-i 

J±(a) = / ^~ yd + t)dt analytic for - oo < Re a < \, (70) 

Jo C 

and (cf. eq. (|21|) ) these integrals admit meromorphic extensions to the whole plane, with 

•11, n f -A 1 (lQg|Cl) (w) (|) , 9 (7U 

simple poles at a = — , of residues ; rH-, n = 1, 2, . . . (71) 

P P 2' 2 (n-1)! 

(the difference J^(a) — J^(a) = 2m2 2a is entire). 

All in all, we finally get two complex conjugate Mellin representations for Z(a): 

-Z(2a) + 2 2ct e* 2 ™ smna [+ e±is °° £' 
Z(cr) = + / r 2a U\ + t)&t; (72) 



2 COS 7TCT 7T Jq £ 2 

and one real principal-value integral representation given by their half-sum, 
= -Z(2.) + 2^2 I% «,r ri ,r 

2C0S7T(J 7T Jo C 2 

(each of the above converges in the full half-plane {Re a < |}). 

Another real form can be obtained with a regular integrand, directly from eq. (|S7l) : 



2 COS 7T(T 7T 



+ + (74) 



however this last integral only converges in the strip {0 < Re a < |}. 
Remarks: 

- as analytical extension formulae, eqs. (JZ2HHI) are precise counterparts of the func- 
tional equation for ((s); they also stand as more explicit forms of Guinand's functional 
relation for Z(a) [T7j, as discussed below (eq.(J7§J)). 

- a similar formula exists for the function £(s,x) of eq.(j3J ([HI], middle of p. 149), 
only requiring Re x > 1 (which precisely avoids the problem raised above by the pole 
of C( s ))j m comparison, the present results correspond to the fixed value x = \, since 
Z(a) = (27r)- 2CT (2cos7ra)- 1 ^(2(T,i) by eq.©; 

As we will elaborate next, analytical properties of Z(a) in the half-plane {Re a < -} 
are made totally straightforward by the Mellin formulae (|72H74j) (while Z(a) is holomor- 
phic in the half-plane {Re a > |}, where its defining series (J58j) converges). Detailed 
results are also recollected in fully reduced form in Sec. 7, Table 1. 
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5.3 Properties of Z{cr) for Re a < 1 

A few leading properties of Z(p) in the half-plane {Re a < 1} emerge more easily as 
special cases from Sec. 4 (although they can be drawn from eq. (l72|) as well): 

- a = \ is a double pole, with the same full polar part as for Z(a), eq. (|^0|) : 

- by specializing eqs. (J5T|) . (J53)l . 

Z (0) = Z{0) = I; Z'(0) = Z'(0) - logH(§) (« 0.811817944). (75) 

Otherwise, a Mellin representation like (|72j) gives a better global view of Z[o) over the 
whole half-plane {Re a < |}. Indeed, its non-elementary part (the integral (|7U|) ) becomes 
regular there, hence can be ignored both for the polar analysis and (thanks to the singer 
factor) for the "trace identities" at integer a: all such information lies then in the first 
term alone, accessible by mere inspection. We thus obtain that: 

- Z{cr) only has simple poles at the negative half -integers a = | — m, with residues 

n m = [ —^[Z(l - 2m) + 2 1 " 2 " 1 ] = - 2 l - 2m )B 2m , m = 1, 2, . . . (76) 

Z7T o7rm 

(hence, only the leading pole a = \ stays double); 

- at the negative integers, the "trace identities" read as 

Z{-m) = L_J_[-Z(-2m) + 2" 2m ] = (-l) m 2" 2m (l - -£ 2m ), m = 0, 1, . . . (77) 

Z o 

(both formulae (|75 j) . (|77j) were fully reduced using eq.((BDJ)). 

- a (cr — > — oo) asymptotic formula follows for Jjjr(er), from the term-by-term substi- 
tution of the Euler product for £(s) into the integrand of eq. (|?2"j) . giving 

Jf(rr) ~ -T(l -2a)J2 A (n)n- 1/2 (\ogn) 2u - 1 , a -> -oo (78) 

n>2 

where as usual, A(n) = logp if n = p r for some prime p, else 0. (An asymptotic formula 
for Zip) itself then follows from eq. (|72*jl and the functional equations for ((s), P(s).) 
Remark: in our notations, Guinand's functional relation for Z(<j) JZj reads as 

2cOS7T<7 2r(2a) COS7TO" 

where Z p (l -2a) ^lim | ^ A(n)n -1/2 (log7i) aT - 1 - / e^ 2 x^Mx} — subject to the 

Riemann Hypothesis jH| — clearly specifies a (real- valued) resummation of the divergent 
series in eq.(f7Sj) (the asymptotic series for —Jp(a)/ r(l — 2a)). Eq.(|79*|l was only asserted 
for < Re a < ~, with no clue as to the analytic structure of either Zip) or Z p (l — 2a) 
elsewhere. The present formulae ([72H74)) are thus resummed versions of eq.(|79|). with a 
definitely more explicit content. 
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5.4 Properties of Z{p) for Re a > - 

As stated before, Z(p) is regular in the half-plane {Re a > -}, where analytical results 
are identities directly obtainable by expanding the logarithm of the functional relation 
()57|) in Taylor series around t — 0. Here we will extract those results from the Mellin 
representation ([72)1. invoking the meromorphic properties of its integral term in the whole 
plane as given by eq. (f7Tj) . 

- for half-integer a = | + m: the residues of the two summands in (f72|) have to cancel 
given that Z{cr) is analytic in the half-plane; this imposes 

(log|C|) (2m+1) (|) = (2m)! [Z(2m+l) + 2 2m+1 ], m = 1,2,... (80) 

( = -2- 2m - 1 (io g r) (2m+1) (i) ), 

which simply amounts to (logH)( 2m+1 )(|) = (itself a consequence of the functional 
equation S(| + i)=H(| — t)); that result further reduces, using eq. lJoTj) . to the identity 

(log|C|) (2m+1) (i) = \{2m)\ (2 2m+1 - l)C(2m+ 1) + \v 2m+1 \E 2m \. (81) 
The case m = is singular, but (logH)'(i) = directly yields 



£ V2/ 2 



log--{lo K n'(i) =Ilog87r + ^ + | (^ 2.68609171). (82) 



- for integer a = m, the pole of the integral is cancelled by the zero of sin na, and the 
following explicit relation results, 

2(-l) m+1 Z(m) - 1 (log |C|) (2m) (|) = Z(2m) - 2 2m , m = 1, 2, . . . . (83) 
(2m — 1)! 

which can also be further reduced with the help of eq. ()61|) . see Table 1. 

Unfortunately, we hardly know anything else about the values (log |C|) < - n - ) (|), n = 
0,1,... (cf. j2H] for C(|))- To supplement the relation (|8*Tj) with ((n) for n odd, we 
can only refer to other formulae for ((2m + 1) (compiled in [S]), and to Euler-Maclaurin 
formulae for Q n '(s) (valid at s = |) with related numerical data pHIHlE]- So, even at v = 0, 
the transcendental values d a Z(0,v) (eq.((T3J)) and Z(m,v) currently remain more elusive 
than at the (exceptional) point v = \ (Sec. 3. 3). Furthermore, we found no reference at 
all to those values (i.e., Z'(0), Zijn)) in the literature. 



5.5 Speculations and generalizations 

The results of Sees. 5. 2-4 for Z((i) are similar to those yielded by the "sectorial" trace 
formula for the analogous spectral zeta function Zx(&) over a compact hyperbolic surface 
X [7|EHI- The present formulae for the Riemann case nevertheless show several distinctive 
features. 

- as announced end of Sec. 2.1, the sequence {t^ 2 } and the analogous spectrum of the 
Laplacian on X have mutually singular features: the former has the parameter values 
/i = |, r = 2 (Z(cr) has its leading pole double, at a = |), whereas the latter more 
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precisely has [1q = 1, r = 1 (Zx{c) has all its poles simple, starting at a = 1), hence this 
spectral analogy for the Riemann zeros holds only partially; 

- in the continuation formulae (JZ2H3J), C( s ) itself reenters as an additive component of 
the shadow zeta function Z(s). This is an altogether different incarnation of ((s) from 
its initial, multiplicative involvement, which remains in the integral term and indirectly 
through the zeros, in the left-hand side Z(o). It is curious to find two such copies of ((s) 
to coexist in one formula, especially with the additive C{2a) represented in its critical 
strip; 

- however, those formulae relative to ((s) are not fully closed as they also invoke the other 
Dirichlet series /3(2cr) (as second additive component in the shadow zeta function Z(2cr)). 
The question then arises whether (3{s) and other zeta functions can be handled on the 
same footing as ((s) (as in [21]), so we now outline a possible extension of eqs. lff^HTIj) . 

We assume that ((s) is a Dirichlet zeta or L-series, having: 

- a single pole, at s — 1 and of order q (typically, q = or 1); 

- the asymptotic property logC(s) = o(s~ N ) for all iV (s — > +oo); 

- a functional equation of the form 

A(t 2 )^|^ = C(| + t), where (84) 

- A{t 2 ) is an entire function of order < 1 in the variable t 2 , and 

- D(t) is an entire function with all of its zeros lying in the half-plane {Re t < 0}. 

Then the Zeta functions Z^ (c) (for the zeros of A) and Z^(er) (for the zeros of D) are 
related by this formula corresponding to eq.(J71J) (we omit the others), 



Z f -(2<r) sinvra f + °° 9<T rC',i x Q 

1 r - f i+t) + -^- T dt. (85) 



s l COS 7TCT" 7T 



2 



Apart from £(s) itself, with eq . ()74p . the next independent example is /3(s). Its func- 
tional equation (fT4"j) has the form (|54*|) with q = (no pole) and A(t 2 ) = H X4 (| + t), 

D(t) = (f)^ + Vr(f + |); its spectrum of trivial zeros (for j3{\ - t)) is {-\ + 2k} (= 



3 7 
2' 2' 



giving as shadow zeta function 



M°) = EH + 2k )' S = 25 [!((! - 2_S )C(^) - • (86) 

fe=i 

Under q = 0, all Mellin representations ([72H74)) coalesce into the single regular form 
. . Z«(2cr) sin7rcr f + °° o„0' , i , , 1s 

and all consequences previously drawn for Zip) have analogs for Zp{u). 

Various such integral representations will naturally add, whenever the initial zeta 
functions combine nicely under multiplication. For instance, eqs. ()74|) and (|87|) add up to: 



2cOS7TCT 71" ' 1 ' % ' ' 







2 



dt. 
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Here, the shadow zeta function purely invokes ({s); on the other hand, under the integral 
sign we now find [^](| + t) so that the new multiplicative zeta function is f3(s)((s), 
also recognized as \ times Z^s), the zeta function of the ring of Gaussian integers Z[i] 
0. Hence eq. lJHHj) becomes 

, , 2 2(J [(1 - 2- 2<T )C(2a) - 11 sinvra f + °° 9rT [Z 4 ' n , In, 

2 COS 7T<7 71 Jo L Z 4 z t — i J 

thus, to isolate C( s )> nere m the additive position, we again had to allow a different zeta 
function elsewhere, this time Z±(s) in the multiplicative position. 

Likewise, by subtracting eq. lfKTJ) from (|71j) instead, we could get the shadow zeta 
function to be /3(s); then the counterpart of eqs. (|ST|) . (jH^j) is a fully explicit identity, 

(log |C|) (2m+1) (|) " (log/3) (2m+1) (|) = ln 2m+1 \E 2m \ + 5 m , log2, m = 0, 1, . . . , (90) 
whereas each of the two left-hand-side terms separately needs ((2m + 1) (or 7 for m = 0). 



6 More about the Hurwitz-type functions 

The purpose of this Section is twofold. First, we analyze the Zeta functions Z{c,v) = 
z 7 2k( T k 2 + v )~ u m ore explicitly over the half-plane {Re a < 1} than in Sec. 4, by exploiting 
the latest special properties of the function Zip-, 0) (with new results even for the case 
v = j). Then, by the same approach, we (briefly) discuss the other Hurwitz-type Zeta 
functions 3 (a, a) and £(s,x), defined through eqs.© and (HJ) respectively. 

6.1 Further properties of Z((7,v) for Re a < 1 

To describe the Hurwitz-type function Z(cr, v) better, we now systematically expand it in 
terms of Z(<r), as 

OO OO / \ 

z(o, v) = v (i +^y°=H p w _7: T) z(°+d '•' (h<v). (si) 

Such an expansion can be formulated around any reference point Vq, but it will be specially 
useful for v = as above. For instance, coupled with eq.(|72j) (say), it can express the 
meromorphic continuation of the general Z(cr, v) to {Re er < |}, while we lack an analog 
of eq. (J72|) itself for any v ^ 0. 

For the polar structure of Z(cr, v ) at a = —m + |, m G N, the series fl9*T|) reduces to 

m Y(- + m-e) 

Z(-m+±+e, v) = — - — ~j~— r Z{-m+£+±+e) v e [+regular part for e -»• 0]; 

2 ' ^ Z\T(\ + m- 1 - e) 2 

(92) 

then, importing the polar structure of Zip) from eq.(|7H]l. we get 

Z(-m + l + e,v) = K — ^ v m e' 2 + TZ m (v) e~ l + 0(1) (e -> 0), (93) 
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just by brute- force polar expansion of the right-hand side of eq.((52jl. Here, the polar part 
of order 2 at every (— m + |) is clearly induced by the only such part of Z{p) (at a = |), 
through the term with i = m in eq . (|92j) ; whereas the residue TZ m (v) is built from all 
residues of Zip) at poles with a > —m + |, as 

(the residues 7?.^ of Zip) at — j + | are known from eq.(J7BJ)). 

Remark: for m — 0, the full polar part (f4*Uj) at cr = |, independent of v, is recovered. 
When a G — N, the series (pTTj) also terminates, as 

m / \ 

Z{— m, v) 

= zZ(p )2(-m + £)v l (meN), (95) 

so that explicit "trace identities" for general v derive from those for v — (eq.(J77J)). (For 
f = j, this result simplifies further, see Table 1.) 

Remark: in view of eqs. ()23j) and ()25j1 . the latter two results now imply a general-n 
formula for the coefficients 0(i_ n )/2(u) in the large-A expansion (fTKJ) of log A(A | {rfe 2 + t>}). 
(Hitherto we had such a formula just at v = 0, not even at t> = |, and knew only the 
other coefficients S(i_ n )/2(u) for any t> and n, by eqs. ([5THo^|l .) 

Our initial emphasis on the special case v = | might now seem misplaced: why didn't 
we operate at once from t> = ? In the first place, we saw the case v — ■? arise more 
readily from the standard product representation of ({s). But mainly, the case v — \ 
also enjoys certain special properties, this time with the values Z'(0) and Z(n) (Sec. 3. 3), 
and since these evolve from transcendental functions of v (Sec. 4), their expansions (|91|) 
around v — are now infinite. So, each case v — and v — | has its own exceptional 
features, the former in the half-plane {Re a < 1}, and the latter for a e N. 



6.2 The Hurwitz-type functions 3(c, o) and £(s, a?) 

The function 3(c, a) as defined by eq.© is a priori more singular than Zip id) (the 
sequence {rfc} itself has r = 2 and /xo = 1, which would require a formalism more elaborate 
than in Sec. 2). Fortunately, 3(o r , a) can also be analyzed directly through its expansion 
around 3(c, 0) = Z(cr), by analogy with eq. (f9"Tjl (see also |S]): 

3(a, a) = (l + f ) ^ = E wl 2ff _ ^ ^ + §4 °* (H<n). 06) 

fc=0 fc £=0 ^ ' 

This formula generates a pole for 3(c, a) now at every half-integer |(1 — n), n G N, 
according to: 

3(|(1 -n) + £ ,a) = ]P ^ r ^ W _~ ^ - n + I) + e) a" +0(e). (97) 
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Differences with eq.()92j) arise due to the factor T{n — 2e)/ T{n — £ — 2e) vanishing whenever 
t > n > 0. Only the polar part at a = ~ remains the same as for Z[p, v) (of order r = 2 
and independent of a, given by eq.()40p): all other poles ~(1 —n) of 3(c, a) are now simple, 
of residues 



1 „ ( n — 1 

r n (a) = — — a 
Aim 

0<2m<n 



+ £ ( J L_ 1 i) W ™ flB " am ' - = 1,2,..- (9f 



(again, lZ m is the residue given by eq. 1)76)1 ). In addition, at a = the e-expansion of 
eq.(jnZJ) captures the finite part too: 

Ct 7 lo 2 ! 71 

ri(o) = Res CT=0 3(o", a) = — -r- ; finite part: FP CT=0 3(er, a) = - H — a. (99) 

47T 8 27T 

As for the function £(s, x), if we express it by eq.© in terms of 3(c, a), then we find 
this combination to be less singular overall: by mere substitution of eq.()97j). £(s, x) shows 
a simple pole at s = 1, of residue — 7r |54"j . and all other possible poles at s = 1 — n, n = 
1, 2, . . . cancel out, resulting in the holomorphy of £(s, x) for all s ^ 1 with the computable 
finite values ("trace identities") 



i+y) 



(27T)"- 1 



n~2m— 1 



71= 1,2,.. 



-7T- 

i n 

0<2m<n 

(100) 

(An alternative evaluation follows from Deninger's continuation formula for £(s,x) (|13j. 
middle of p.149), as f (1 -n,\ + y) = {2>K) l - n [{y + i)^ 1 + (y - i)™" 1 + 2™- 1 5 n (± + ±y)], 
whose agreement with eq. (jlOO)) can be verified.) 

As for special values: first, <9 s £(s, x) s=0 is expressible as well, in terms of £(x) [T31 I34"]: 

d£(s,x) s=0 = log 2 1 / 2 (2tt) 2 - log ~(x) (101) 

<^=> -<9 s [^(x - p)- s ] s=0 = logH(x) + i(log27r) x - ilog47r (102) 
p 

(the equivalence of the two forms follows from eqs.(|4f and (jlOO)) for n = 1, i.e., £(0,x) = 
|(x + 3)). Now, the exponentiated left-hand side of eq. ()102|) precisely defines the zeta- 
regularized product A zr (x) built upon the sequence {p} of Riemann zeros, while the right- 
hand side mainly involves S(x) of eq.(|28)l. So, eq. p()2jl is converting a zeta-regularized 
product (Azir(^)) to Hadamard product form. As an aside, we now verify that such a 
conversion formula is entirely fixed by universal rules for (complex) admissible sequences, 
specialized here to r = 1 (as in [HUGH]) and p = 1 — since the Zeta functions £(s, x) 
have just a simple pole at s = 1. Those rules yield these two prescriptions: log A zr (^) = 
logH(x) — (ax+/3), and the large- x expansion of log A zr (^) shall only retain canonical (or 
standard) terms, namely: c^x^ for 1 > /j ^ N, cix(logx — 1), cologx. Those conditions 
together fix [a, (3) uniquely, and here, eqs. (JHSHSU) for logS(x) as input precisely lead to 
eq. ()102|) as output. 

Likewise, the special values x), n = 1, 2, . . . are expressible in terms of [the higher 
Laurent coefficients of] C( x ), e -g-, by applying residue calculus to Deninger's continuation 
formula ([El, p. 149). 

Thus, a fair degree of structural parallelism finally shows up between the two Hurwitz- 
like families Z(s,v) and x) 
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^)=E(r, 2 + l)- CT [v 

k=l 



Z{a) = Y,T k - 2 ° [v = 0] 
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Z'(0) = log[2 11 /4 7r i/2 r( i ) -i| C( i; 
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Table 1: Analytical results for £-Zeta functions of the Riemann zeros. Notations: see 
eqs.dH),©, (|TT|) [for /3(s)], (}4T)j) [for 7£]; ^ m = J2 P P~ m 'i m stands for any positive integer, 
and e — > 0. As an indexing tool, the superscripts refer to the relevant equation numbers 
in the main text. [In the very last formula (for Z(m)), one can still use eq.(|HJ) for ((2m)]. 



7 Recapitulation of main results 

In way of conclusion, Table 1 collates the analytical results found for the two £-Zeta 
functions Zip) (= Zip-, v = |)) and Zip) (= Zip, 0)). Furthermore, corresponding results 
for the general Z(<y,v) were derived in Sec. 4 (for Re a > 0) and 6.1 (for Re a < 1), and 
partly extended to the functions 3(a,a) and £(s,x) in Sec. 6. 2. 

The other novel results we have developed here concern Zip) in the half-plane {Re a < 
|}: the analytical continuation formulae (J72HHJ); and the a —* — oo asymptotic formula 
(J78|) as corollary. We also came across two elementary (but unfamiliar to us) formulae 
concerning ((s) itself: eqs.()81| ) .([90| ) . 

Appendix B gives information on some numerical aspects based on our use of the 
100,000 first Riemann zeros (made freely available on the Web by A.M. Odlyzko [3U], to 
whom we express our gratitude). We also wish to thank C. Deninger, J. P. Keating, P. 
Lebceuf, V. Maillot, C. Soule, and the Referee, for helpful references and comments. 
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Appendix A: Meromorphic Mellin transforms 

We briefly recall the meromorphic continuation argument for a Mellin transform like 

eq.((Tni), 1(a) == J +o ° L(z)z~ a ~ 1 dz, assuming the function L(z) to be regular on IR + (for 
simplicity), with 

oo 

L(z) = O(z wo ) (z -> + ), and L(z) ~ ^« log ,2 + a^Jz^ (z -> +00) (103) 



n=0 



as in eq. (|TK|) (asymptotic estimates are repeatedly differentiable); and crucially, /io < vo- 

Sequential directed integrations by parts can be used (see j2UE3El for details). 
Step 1: 1(a) converges for /i < Re z < u , and in that strip, 

r+oo 7 IM>-<7 

1(a) = [L(z)z-^\ dz. (104) 

Jo o - /i 

If a M0 = 0, this suffices: the new integral actually converges for < Re z < is (thanks 
to z fl0+1 [L(z)z^ tJ ' ]' = 0(z t11 log z) for z — > 00), hence /(a) is manifestly meromorphic in 
that wider strip, with a simple pole at a = // of residue 



[L(^)^ Q ] , d^ = a Ai0 ; (105) 

furthermore, in the complementary strip {/ii < Re z < fio}, backward integration by 
parts now yields 

r+00 

1(a) = / [L(z)-a^)]z~°- l o\z. (106) 







Then the whole argument can be restarted from here, to extend 1(a) further (across 
{Re a = /ii}), and so on: the case r = 1 thus gets settled. 

Step 2: if a w 7^ 0, one more integration by parts upon eq. (|104|) yields 

r+OO yW-o 

1(a) = / [z[L(z)z-^]'}'- -dz (jii < Re z < Uq). (107) 

Jo O - M 1 

All previous arguments then carry over, but the pole is now double, with 

r+00 

principal polar cefficient: / [z[L(z)z^ fl0 ]']' dz = a^ , (108) 

Jo 

and residue — J +oo [z[L(z)z~^ }'}' log zdz (from the residue calculus); integrations by parts 
(backwards, and split) reduce the latter to 

rl r+00 1 

residue: / [L(z)z~^]' dz + / {z[L(z)z~^ }' - a^}- dz = a m . (109) 
Jo Ji z 

The last two formulae thus generate eq. (|20j) for the leading double pole, and so on for 
r = 2. 
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(More generally, if the factor of z^ n in the expansion ()103j) is a polynomial of degree 
p n in log z, then /i n becomes a pole of order (p n + 1) for 1(a) [2*T].) 

As for meromorphic continuation in the other direction (across {Re a = vo})-, it works 
likewise if L(z) admits a z — > expansion: the previous arguments apply upon exchanging 
the bounds and +oo under cr i— > —a. E.g., in the regular case, L(z) expands in an entire 
series at z = 0, and step 1 suffices (as in the main text, where v = 1). 



Appendix B: Numerical aspects 

We complete our analytical study by describing some very heuristic numerical work with 
Z(a) and Z(cr), mostly in the range {a > 0}, and focusing on the simpler case of Z(c). 
(The same ideas apply for any generalized Zeta function Z(cr, v) and for complex cr, but 
the formulae get more involved.) Here, the Riemann Hypothesis is de facto implied 
throughout (there being no numerical counter-example). 

Numerically, Z(o) looks almost indistinguishable from Z(o) for a > 0, because T\? ^> 
7 (VA;) (an empirical fact; already, r x 2 « 199.790455). If we expand Z(a) in terms of 
Z(o) according to eq. (j9~Tj) . and make the roughest approximations, we get that Z(o) « 
Zip) [1 — "/(4r 1 2 )]: i.e., the very first correction term is only of relative size ~ cr/800. As 
other related numerical observations: 



< Z(o) - Z{a) < 0.0003 for all real a > 0; (110) 

A = 4[Z'(0) - Z'(0)] w -0.0231003495 vs B = -Z(l) w -0.0230957090 (111) 

(not only is A small, but moreover, A = 41ogH(|) by eq. (|73|) and B = [logS]'(0) by 
eq. (|3ip . hence \B — A\ < 5 x 10 -6 reflects how little the function logH(s) deviates from 
the parabolic shape As(l — s) over the interval [0, 1]). 

We now focus on the numerical evaluation of Z(o) itself. The defining series ()58|) 
converges more and more poorly as a —* | (with divergence setting in at o = |). We 
then replace a far tail of that series (k > K) by an integral according to the integrated- 
density estimate (J29|) . and formally obtain a kind of Euler-Maclaurin formula, 

K-l 



Z(o)= lim S K (a), W^Vr^ + ^ + M, (112) 

1 

i t i—2(7 r _ -i 



(similar formulae can be written for Z'{<r), Z(cr,v), finite parts at cr = |, etc.). 

The approximate remainder term Rk((t) balances the dominant trend of the partial 
sums Yl T k~ 2a - It thereby accelerates the convergence of the partial sums in eq. (|112|) 

k<K 

for a > while for cr < | it counters their divergent trend, so that Sk(c) converges (as 
K — > +oo) when a > ([Hj, p. 116, last line). The next obstruction to convergence arises 
at (j = but is of another type: Sk{&) displays erratic fluctuations in K (roughly of 
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oo 


oo 


a 


r- , / \ v — / O 1 \ ^ r In 

Z(*)= E(^ 2 + |)- CT b = |] 


Zip) = E r k - 2(T [v = 0] 

k=l 


-1 


-0 0625* 

V 7 • V J V7 .7 


-0.28125* 


-3/4 


1.69388 


0.54319 


-1/4 


0.800805 


0.785321 





0.875* 


0.875* 


derivative at 


U.oUbU4zy 


O.ollol li) 


+ 1/4 


1.548829 


1.549060 


finite part at 1 /2 


0.251546 


0.251637 


+3/4 


0.247730 


0.247760 


+ 1 


0.0230957 


0.0231050 


+3/2 


0.0007287 


0.0007295 


+2 


0.0000371 


0.0000372 



Table 2: Numerical values for £-Zeta functions of the Riemann zeros (*: exact values). 
Implied precision is expected to hold, but not guaranteed. 



the order r^ _2cr (loglogT^;) 1//2 , according to [SHI, eq.(2.5.7)), and those numerically blow 
up indeed (as K — > +oo) when o~ < 0. Further convergence now requires to perform a 
damping of those fluctuations (as argued previously for "chaotic" spectra [3]). Here, a 

dcf K 

Cesaro averaging (defined by (S)k = K -1 E^') ^PP^ars to work well initially (results 

i 

can be verified at o = 0), but not very far down: already at o = —0.25, the fluctuations of 
(S)k(o') itself retain a standard deviation > 10~ 3 up to K rs 10 5 . So, instead of pursuing 
ever more severe (and unproven, after all) numerical regularizations as a decreases below 
i, we advocate the switch to the continuation formulae (J72HZ3) for numerical work as 
well. Thus, we first tested eq.fT73^) against eq. (|112|) for Z(+\), then used it to evaluate 
Z(—j), plus eq. lJDTj) with v = ~ (3 terms sufficed) to obtain Z{— ~). 

Table 2 gives a summary of the numerical results we obtained. (We found no earlier 
analogs, except for the other special sums ^ n in J2H]-) 
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